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Abstract. The 1-Jettiness (ri) event shape for Deep Inelastic Scattering (DIS), allows for a 
quantitative and global description of the pattern of QGD radiation for single jet (J) production 
in electron-nucleus (Na) collisions e~ + Na e~ + J + X. It allows for precision studies 
of QGD and is a sensitive probe of nuclear structure and dynamics. The large transverse 
momentum {Pjj,) of the final state jet J, characterizes the hard scale in the problem. The 
region of phase space where ri <C Pj ^, corresponds to configurations where energetic radiation 
{E ~ Pjt) is only along either the single jet direction or the beam direction with only soft 
radiation (i? ~ ti <C Pjt) in between. Thus, the restriction n <C Pjj- corresponds to a veto 
on additional jets and leads to large Sudakov logarithms of rijPj^ that must be resummed. 
Based on a factorization framework, derived using the Soft Gollinear Effective Theory (SGET), 
we provide resummation results at the NNLL level of accuracy and match them onto the NLO 
result in fixed order perturbation theory, appropriate in the ri ~ Pj^ region where additional 
jets and hard radiation are allowed. The ri-distribution depends on the jet algorithm used to 
find the leading jet in the region n ~ Pj ^, unlike the resummation region where this dependence 
is power suppressed in TxjPjrj. 1. We give results for the entire n spectrum, with a smooth 
matching between the resummation region and fixed-order region, where we make use of the 
anti-kt jet algorithm. The 1-Jettiness event shape can be a powerful probe of nuclear and QGD 
dynamics at future electron-ion colliders and by analyzing existing HERA data. 


Event shapes probe QCD dynamics, providing a quantitative global measure to characterize 
QCD radiation in the final state and allowing good analytic control over the corresponding 
calculations. Various event shapes for Deep Inelastic Scattering (DIS) processes were first studied 
in Refs. [U |2l 131 H] . Thrust m and Broadening [3] distributions were studied at the next-to- 
leading logarithmic (NLL) level of accuracy and matched at 0(as} to fixed-order results. A 
numerical comparison was also done against ©(a^) results [5l |6]. Thrust distributions were 
measured at HERA by the HI [71 [8] and ZEUS P collaborations. 

In Refs, [m [HI |T3l IHj , single jet (J) production in the DIS process 

e- +Na^ J + X, (1) 

^ ^Speaker at PAVI 14: From Parity Violation to Hadron Structure, Skaneateles, NY. 




Figure 1. Typical configuration Figure 2. Typical configuration 

for Ti <C Pjj ,. for Ti ~ -Pjj,. 


where Na denotes a nucleus with atomic weight A, was first studied using the 1-Jettiness event 
shape (ri), a specific application of the A^-Jettiness event shape [l5l [H] first introduced to 
study exclusive A^-jet production at the LHC. A factorization and resummation framework was 
derived mm for the observable 
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in the limit ri <C Pjj,, where Pjj, and y denote the transverse momentum and the rapidity of 
the jet (J). The 1-jettiness global event shape n is defined as 
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where the sum is over all final state particles (except the final state lepton) with momenta 
Pk- The light-like four-vectors qa and qj denote reference vectors along the nuclear beam and 
final state jet directions respectively. In general, an external jet algorithm is used to determine 
leading jet and the light-like vector qj is aligned with it. The constants Qa and Qj are of the 
order of the hard scale and their choices are not unique; different choices correspond to different 
definitions of ri. The 1-Jettiness algorithm associates all final state particles either with the 
beam region or with the jet region according to the minimization condition in Eq. Q. The 
momentum of the final state jet Pj, defined in the 1-Jettiness framework, is then given by the 
sum of the momenta of all particles in the jet region 
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Note that the external jet algorithm is only used to determine the light-like reference vector qj 
and that the 1-Jettiness jet momentum Pj is in general distinct from that of the leading jet 
found by the external algorithm, as explained in detail in Ref. m- 

The limit ti <C Pjj, corresponds to configurations that typically look like that shown in 
Fig. 0 any energetic radiation (E ~ Pjj,) in the final state is closely aligned either along the 
jet direction or along the beam direction. At wide angles from these directions, the restriction 
Ti <C Pjj, only allows for soft radiation (E ~ ri <C Pjj-)- In effect, the restriction ri <C Ejj, 
acts as a veto on additional jets or hard radiation at wide angles from the beam or leading 
jet directions. This restriction on final state radiation gives rise to large Sudakov logarithms 
of the form a” In^™ (ri/Pj^) for m < n, that require resummation. Since the dynamics in the 










Ti <C Pjj, region is dominated by radiation collinear with either the jet or beam directions and 
soft radiation in all directions, the Soft Collinear Effective Theory (SCET) [T71 [THlITni 
is the appropriate effective theory to derive a factorization and resummation framework. This 
factorization and resummation framework was first developed in Refs. mm and has the 
schematic form 

da,esum = ~ H0B0J0S, (5) 

dyaPj^dTi 

where H is the hard function, B is the beam function mm that describes the dynamics of the 
initial state PDF and the perturbative initial state radiation collinear with the beam direction, 
J is the jet function describing the dynamics of the collinear radiation in the final state jet, and 
S is the soft function describing the dynamics of soft radiation (E ~ ti) throughout the event. 
The beam function is matched onto the standard PDF B ~ X (D /, where I is perturbatively 
calculable and describes perturbative collinear radiation along the beam direction. The hard, 
jet, beam, and soft functions have renormalization scales with respective scalings 


d-H ~ Pjt^ 
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that minimize any large logarithms in the respective functions. All objects in Eqs. are 
evaluated at a common scale [i after using renormalization group equations in SCET to evolve 
them from their natural scalings in Eq. (§, thereby resuming large logarithms. We refer the 
reader to Refs. mm for a detailed version of the factorization formula, including the derivation 
and field-theoretic definitions of the hard, jet, beam, and soft functions. 

Numerical results with a resummation of Sudakov logarithms at the next-to-next-leading 
logarithmic (NNLL) level of accuracy were first presented in Ref. [12]. These results also included 
a wide range of nuclear targets: proton. Carbon, Calcium, Iron, Gold, and Uranium. Shortly 
thereafter, NNLL resummation results for a proton target were presented in Ref. jUj and they 
also introduced two new definitions of 1-jettiness, corresponding to different choices of the jet 
reference vector qj used in the definition of the 1-jettiness event shape. 

The situation is quite different in the region 


n ~ Pjt, (7) 

for which the configurations typically look like that shown in Fig. The large value of ti 
corresponds to easing the veto on hard radiation at wide angles from the beam and leading 
jet directions. The Sudakov logarithms of ti/Pj^ are now small and the use of standard 
perturbation theory is appropriate. 

In order to obtain the full ri-spectrum, one must smoothly match the resummation (ti <C 
Pjj,) and fixed-order (ti ~ Pjj,) regions. This was done in Ref. [2], where the resummation 
region was smoothly matched onto the fixed-order region, at the NNLL -|- NLO(~ a^). More 
recently, results at NNLL-I-NLO were also obtained in Ref. |25| for a different definition of 1- 
Jettiness, which was shown to be equivalent to the DIS thrust [I] event shape and does not use 
a jet algorithm in its implementation. Schematically, the differential cross-section for the full 
spectrum can be written as 

da = [dcJresum “ daf£^^] + da^^ . (8) 


Here dciresum denotes the resummed cross section computed in the region t\ <C Tj^. The 
d^fiSnra is this resummed cross section expanded to fixed-order perturbation theory and is given 
by setting all scales in the factorization formula equal to each other 

d<yfe?nra = daresum{^^ = (9) 




thereby turning off resummation and leaving only the contributions of fixed-order SCET matrix 
elements. The is the full cross section at the same order in perturbation theory. The 

differs from by terms that are non-singular in the limit ti —)• 0. In the resummation region 

Ti <C Pjj,, da is dominated by duresum due to a cancellation between and da^^, up to 

suppressed non-singular terms. Similarly, in the fixed-order region ri ~ da is dominated by 
da^^ due to a cancellation between duresum and up to terms suppressed in perturbation 

theory. Furthermore, in order to smoothly match the ri-spectrum in the resummation and fixed- 
order regions, one must make use of profile functions [26l [271 EH] so that the scales hh, I^b, l^j, 
and fis appearing in ddresum smoothly converge to the single scale /r ~ Pj^, that appears in 
da^^. This is essential for important cancellations to occur between the various terms Eq. ([^. 


The full Ti-spectrum has three distinct regions 
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that must be smoothly connected by matching and the use of profile functions. In the region 
Ti ~ ^QCD the soft radiation (E ~ ri) becomes non-perturbative and must be modeled. The 
model must smoothly converge to the perturbative soft function in the region n !§> Aqcd- This 
is accomplished by writing the soft function as a convolution |26[ 129] between a model function 
and the partonic soft function. 

Note that in the resummation region ri ^ the dependence on the jet algorithm used 
to determine qj in Eqs. Q and Q power suppressed [H] in n/Pjj,. This can be understood 
by noting that ri <C Pjj, corresponds to the configuration in Fig. with a narrow collimnated 
jet of energetic radiation well-separated from the beam direction with only soft radiation in 
between. For such configurations, different jet algorithms will find the same direction qj for the 
leading jet momentum, up to power corrections in ri/Pj^ corresponding to differences in how 
soft radiation is grouped into the leading jet. On the other hand, for ri ~ Pj^ as in Fig. 
different jet algorithms can give rise to different directions qj for the leading jet momentum 
corresponding to differences in how hard radiation at wide angles is grouped into the leading 
jet. 

Thus, in the schematic formula in Eq. ([^, diTresum and are independent of the jet 

algorithm of used to find the leading jet, up to power corrections. This allowed the works 
in Refs. na Ea ED to provide NNLL results in the resummation region ri <C Pjj, without 
reference to an explicit jet algorithm. On the other hand, in the region ti ~ Pj^, an explicit 
jet algorithm must be used and the reference vector qj will strongly depend on the algorithm 
used. In particular, the computation of da^^ in Eq. ^ requires an explicit jet algorithm. In 
our work we make use of the anti-kt m jet algorithm, although our numerical code is flexible 
enough to use other algorithms. The da^^ has the schematic form 

da^^ ~ IdPS E„([E5])|M|2®/, (11) 

where dPS is final state phase space measure, \Ai'^ is the UV renormalized amplitude squared 
for the partonic process, / denotes the initial state PDF, and Emeas. is the measurement function 
that imposes restrictions on the final state. In particular, for the observable in Eq.Q it restricts 
the final state jet to have a transverse momentum and rapidity of Pjj, and y respectively and 
the final state radiation to have the value ri for the 1-jettiness event shape.These final state 
restrictions along with the anti-kt jet algorithm are implemented numerically using Vegas m- 
For each phase space point, a jet algorithm is implemented to cluster final state particles and 
find the leading jet. The transverse momentum {Kj^) and rapidity {yx) of the leading jet are 


then used to construct the light-like jet reference vector qj = {Kj^ coshyK, Kj^,K sinhyi^:). 
A set of values ri, Pjj,,y is returned for each phase space point. Numerical integrations are then 
performed by restricting the phase space to be within specified bin sizes around specified values 
for ri,Pjy, and y. 

The partonic channels LO are 


e +qi^ e + qi, 

e~ +qi^ e~ + qi (12) 

where the index i runs over the quark and antiquark flavors. The NLO contribution has three 
types of partonic channels with the real emission of an extra parton in the final state 


e~ +qi^ e~ +qi+ g, 
e~ + qi ^ e~ + qi + g 

e~ + g ^ e~ + qi +qi, (13) 


and virtual corrections to the leading order channels in Eq.(12). Infrared (IR) singularities arise 
in these NLO calculations from the real emission of an extra parton in the final state as well 
as from the virtual corrections to the leading order process. In order to numerically evaluate 
d(j^^, it becomes necessary to analytically isolate these IR singularities. We use dimensional 
regularization, working in d = 4 — 2e dimensions, to isolate the IR divergences as poles in e. 
We also implement the sector decomposition technique [32l ESI EH ES] in order to break up the 
phase space into sectors where only single parton or a single pair of partons becomes unresolved, 
corresponding to the soft and collinear IR divergences. This facilitates the isolation of IR poles 
after which numerical integration techniques can be used straightforwardly. We refer the reader 
to Ref. [iH for more details. We now give numerical results for the full 1-Jettiness spectrum for 



Figure 3. ri-distributions for 
the full NLO (solid red line) and 
the magnitude of the resummed 
SCET result expanded to NLO for 
kinematics given in the text. 



T|(GeV) 

Figure 4. The NLL (widest blue 
band), NLL’+NLO (red band), and 
the NNLL+NLO (green band) ti- 
distributions. 


the case of a proton target. We work at an electron-nucleus center of mass energy of 90 GeV and 
integrate over the jet transverse momentum and rapidity over the ranges [P^, T’j^®*^]=[20 GeV, 
30 GeV] and \y\ < 2.5 respectively. In Fig. we show the perturbative results for /dri 
(solid red curve) and \daf^-^^/dTi\ (dashed blade curve). We see that in the limit n —>■ 0, the 
da^^ converges to as expected since in this limit da^^ is dominated by terms singular 

in this limit. This is one of several cross-checks performed [TT] on the NLO calculation. 











We also note that around ri ~ 5 GeV, becomes negative. This corresponds to the 

region where the non-singular terms, not contained in become important. A smooth 

matching of the resummation and fixed-order regions, as in Eq. Q is required to properly 
describe the spectrum over its full range. The result of this matching which also incorporates a 
non-perturbative soft function in the region ri ~ Aqcd is shown in Fig. The various bands 
in Fig. [^correspond to the Ti-distribution at the NLL (widest blue band), the NLL’-I-NLO (red 
band), and the NNLL-I-NLO (green band) levels of accuracy. The width of the bands indicate 
the perturbative uncertainty obtained via scale variation. We see that resummation tames the 
singular behavior in the ri —)■ 0 limit and that the perturbative uncertainty reduces as we go to 
higher orders in resummation. 
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